Antifreeze proteins (AFPs) are synthesized by various organisms to enable their cells to survive subzero environment. These proteins bind to small ice crystals and inhibit their growth, which if left uncontrolled would be fatal to cells. The crystal structures of a number of AFPs have been determined; however, crystallographic analysis of AFP–ice complex is nearly impossible. Molecular modeling studies of AFPs’ interaction with ice surface is therefore invaluable. Early models of AFP–ice interaction suggested H-bond as the primary driving force behind such interaction. Recent experimental evidence, however, suggested that hydrophobic interactions could be the main contributor to AFP–ice association. All computational studies published to date were carried out to verify the H-bond model, and no works attempting to verify the hydrophobic interaction model have been published. In this work, we Monte Carlo–minimized complexes of several AFPs with ice taking into account nonbonded interactions, H-bonds, and the hydration potential for proteins. Parameters of the hydration potential for ice were developed with the assumption that the free energy of the water–ice association should be close to zero at equilibrium melting temperature. Our calculations demonstrate that desolvation of hydrophobic groups in the AFPs upon their binding to the grooves at the ice surface is indeed the major stabilizing contributor to the free energy of AFP–ice binding. This study is consistent with available structural and mutation data on AFPs. In particular, it explains the paradoxical finding that substitution of Thr residues with Val does not affect the potency of winter flounder AFP whereas substitution with Ser abolished its antifreeze activity.
展开▼
